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Abstract 

We present an experimental characterization of the statistics of multiple photon pairs produced 
by spontaneous parametric down-conversion realized in a nonlinear medium pumped by high-energy 
ultrashort pulses from a regenerative amplifier. The photon number resolved measurement has been 
implemented with the help of a fiber loop detector. We introduce an effective theoretical description 
of the observed statistics based on parameters that can be assigned direct physical interpretation. 
These parameters, determined for our source from the collected experimental data, characterize the 
usefulness of down- conversion sources in multiphoton interference schemes that underlie protocols 
for quantum information processing and communication. 
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FIG. 1: (Color online) A schematic of a spontaneous parametric down- conversion source. The 
nonlinear crystal X is pumped with a laser beam p. Generated photons are highly correlated and 
useful modes a and b are typically selected by narrow spatial and frequency filters /. 



I. INTRODUCTION 

Spontaneous parametric down-conversion (SPDC) is the basic source of non-classical 

n a 

light in experimental quantum optics p|, testing foundations of the quantum theory [2|, and 
implementing protocols for quantum information information processing and communication 
The essential feature of SPDC is the guarantee that the photons are always produced 
in pairs, and suitable arrangements allow one to generate various types of classical and 
quantum correlations within those pairs. 

The physics of SPDC depends strongly on optical properties of nonlinear media in which 
the process is realized. This leads to an interplay between different characteristics of the 
source and usually imposes trade-offs on its performance. For example, many experiments re- 
quire photon pairs to be prepared in well-defined single spatio-temporal modes. In contrast, 
photons generated in typical media diverge into large solid angles and are often correlated in 
space and time, as shown schematically in Fig. [TJ Specific modes can be selected afterwards 
by coupling the output light into single-mode fibers and inserting narrowband spectral fil- 
ters. However, it is usually not guaranteed that both the photons in a pair will always have 
the matching modal characteristics, and in many cases only one of the twin photons will get 
coupled in This effect, which can be modelled as a loss mechanism for the produced light, 
destroys perfect correlations in the numbers of twin photons. These losses come in addition 
to imperfect detection, and can be described jointly using overall efficiency parameters. 

The effects of losses become more critical when the SPDC source is pumped with powers 
so high that it is no longer possible to neglect the contribution ot events when nrultip.e pairs 
have been simultaneously produced |5|. Such a regime is necessary to carry out multiphoton 
interference experiments, it can be also approached when increasing the production rate of 
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photon pairs. One is then usually interested in postselecting through photocounting the 
down-conversion term with a fixed number of photon pairs and observing its particular 
quantum statistical features [^. In the presence of losses the same number of photocounts 
can be generated by higher-order terms when some of the photons escape detection. However, 
the statistical properties of such events can be completely different, thus masking the features 
of interest. Although some quantum properties may persist even in this regime, with a 
notable example of polarization entanglement their extraction and utilization becomes 
correspondingly more difficult. 

The present paper is an experimental study of multiphoton events in spontaneous para- 
metric down-conversion with particular attention paid to the effects of filtering and losses. 
The multiple-pair regime is achieved by pumping the nonlinear crystal by the frequency- 
doubled output of a 300 kHz titanium-sapphire regenerative amplifier system. The kilohertz 



repetition rate has a 
of the loop detector 



lowed us to count the number of the photons at the output with the help 

n 

|. Using a simplified theoretical description of the SPDC source we in- 
troduce effective parameters that characterize its performance in multiphoton experiments. 
The obtained results illustrate trade-offs involved in experiments with multiple photon pairs 
and enable one to select the optimal operation regime for specific applications. 

This paper is organized as follows. First we describe a theoretical model for SPDC 
statistics in Sec. [Tll Sec. [1111 introduces effective parameters to characterize SPDC sources. 
The experimental setup and measurement results are presented in Sec. HVl Finally, Sec. IVl 
concludes the paper. 



II. SPDC STATISTICS 



We will start with a simple illustration of the effects of higher-order terms in SPDC. 
Suppose for simplicity that the source produces a two-mode squeezed state which can be 
written in the perturbative expansion as |0,0) + r|l, 1) + r^|2,2) + . . ., where r measures 
squeezing and is assumed to be real. For two-photon experiments, the relevant term is |1, 1) 
and the contribution of the higher photon number terms can be neglected as long as r <^ 1. 
This enables postselecting the two-photon term and observing associated quantum effects, 
such as Hong-Ou-Mandel interference. Suppose now that each of the modes is subject to 
losses characterized by l — rj, where rj is the overall efficiency. Losses may transform the term 
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|2,2) into |2,0) or |0,2), whose presence will lower the visibility of the Hong-Ou-Mandel 
interference. The two-photon term now occurs with the probability T^^r^, while the four- 
photon term effectively produces one of the states |2,0) or |0,2) with the total probability 
equal to 2(1 — rjYrfr^. This constitutes a fraction of 2(1 — riYr"^ of the events that come 
from single pairs produced by the source. This fraction can easily become comparable with 
one, especially when the losses are large. 

Let us now develop a general model of photon statistics produced by an SPDC source. In 
the limit of a classical undepleted pump the output field is described by a pure multimode 
squeezed state. By a suitable choice of spatio-ternporal modes, called characteristic modes, 
such a state can be brought to the normal form [9|| in which modes are squeezed pairwise. 
Denoting the annihilation operators of the characteristic modes by and S^, the non- 
vanishing second-order moments can be written as: 

{alai) = {blk) = y(cosh2rfc - 1) 

(afcM = (aISJ)* = ^sinh2rfc (1) 

where is the squeezing parameter for the kth pair of modes. Because the state of light 
produced in SPDC is gaussian, these equations, combined with the fact that first-order 
moments vanish (a^) = (bk) = define fully quantum statistical properties of the output 
field. 

Let us first consider the case when the spatial and spectral filters placed after the source 
select effectively single field modes. The annihilation operators a and b are given by linear 
combinations of the characteristic modes: 

k k 

where tk and t'^ describe amplitude transmissivities of the filters for the characteristic modes 
and ItfcP = J2k l^fcP ~ Because the complete multimode state is gaussian, the reduced 
state of the modes d and b is also gaussian, and it is fully characterized by the average 
numbers of photons n = {d^a) and n' = (b^b)^ and the moment S = {db). These quantities 
can be written in terms of the multimode moments given in Eq. ([T]), but we will not need 
here explicit expressions. It will be convenient to describe the quantum state of the modes 
d and b with the help of the Wigner function: 



where a = {a, a*, (3, ^^^^ C is the correlation matrix composed of symmetrically ordered 
second order moments: 
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(4) 



Given the Wigner function of the reduced state for the modes a and h in the Gaussian 
form, the calculation of the joint count statistics pn,m is straightforward. It will be useful to 
introduce an operator representing the generating function of the joint count statistics 

S(x, y) = ^ m) {n, m\ = x"'^"'y^^^ (5) 



n,m 



whose expectation value over the quantum state expanded into the power series yields the 
joint count statistics Pnm'- 

1 



(6) 



x=y=0 



^"'"^ n\m\ dx^dy 

Because the operator S(x, y) is formally equal, up to a normalization constant, to product of 
density matrices describing thermal states of modes d and b with average photon numbers 
x/{l — x) and y/(l — y) respectively, the corresponding Wigner function has a Gaussian 
form: 



(7) 



71^(1 + x){l + y) 

X exp ^-2|q;|^ TT~) ("^I'^l 
Using the above expression it is easy to evaluate the generating function for the joint 
count statistics by integrating the product of the respective Wigner functions: 



1 + 2/ 



{E{x,y))^7r' / d^ad^(3W^{a,(3)Wp{a,(3) 
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+ 1 - N{r]x + 1 - r]){r]'y + 1 - 77') 
where we introduced the following three parameters: 

-{\S\'-nn') 
' - fin') 
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These three parameters have a transparent physical interpretation in the regime when 
|(a6)p > {a^d){b^b). Then, it is easy to check that the generating function given in Eq. ((HI) 
describes the count statistics of a plain two-mode squeezed state whose two modes have been 
sent through lossy channels with transmissivities 7] and r/', and the average number of photons 
produced in each of the modes was equal to A^. Note that the inequality \ {ab)\'^ > {a)a)(})^h) 
implies nonclassical correlations between the modes a and h. If the opposite condition is 
satisfied, the state can be represented as a statistical mixture of coherent states in modes a 
and h with correlated amplitudes. In this case, it is not possible to carry out an absolute 
measurement of losses. 



III. PARAMETERS 



In a realistic situation, the spectral filters employed in the setup are never sufficiently 
narrowband to ensure completely coherent filtering. Therefore a sum of counts originat- 
ing from multiple modes will be observed. We will model this effect by assuming that the 
detected light is composed of a certain number of modes with identical quantum sta- 
tistical properties described in the preceding section. The generating function E{x,y) for 
count statistics is therefore given by the A^-fold product of the expectation value {E{x,y)) 
calculated in Eq. (IHl): 

Eix,y) = {Eix,y))^ (10) 

The parameter which we will call the equivalent number of modes, can be read out 
from the variance of the count statistics in one of the arms characterized by the generating 
function E{x,y). For a single mode source the variance is equal to that of a thermal state 
with (An)^ = {n) + (n)^. It is easily seen that for M. equally populated modes the variance 
becomes reduced to the value (An)^ = {n) + {n)'^ /M.. Solving this relation for Jvl leads us 
to a measurable parameter that will help us to characterize the effective number of detected 
modes: 

M = (11) 
{Any -in) ^''^ 

Note that Ai is closely related to the inverse of the Mandel parameter [lOj. 

The second parameter we will use to characterize the SPDC source measures the overall 
losses experienced by the produced photons. Let us first note that in the perturbative 
regime, when all the squeezing parameters -C 1, each one of the quantities n, n', and 15*1^ 
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appearing in Eq. ([91) is proportional to the pump intensity. In this regime the efficiencies 
can be approximated by the ratios rj ~ \S\'^/n' and rj' ^ \S\^/n that are independent of the 
pump intensity. On the other hand, the average photon numbers {n) and {n') in both the 
arms are linear in the pump power. 

n 

Following early works on squeezing [11||, we will introduce here a parameter that quantifies 
the subpoissonian character of correlations between the counts n and n' in the two arms of 
the setup. First we define a stochastic variable: 



n n' \ I 1 

¥)~W)] / W)^^)' ^^^^ 

This definition takes into account the possibility of different losses in the two arms through a 
suitable normalization of the count numbers. The subpoissonian character of the correlations 
can be tested by measuring the average (5^). The semiclassical theory predicts (5^) > 1, 
while for two beams with count statistics characterized by the generating function y))^ 
one obtains: 

In the case of equal efficiencies rj = rj' this expression reduces simply to (5^) = \ — rj. We 
will use the last relation as a method to measure the average overall efficiency of detecting 
the state produced by the SPDC source. 

As discussed at the beginning of Sec. [Tll the initial quantum state of light used for many 
experiments should ideally be in a state |1,1) or |2,2). In the perfect case of two-mode 
squeezing and negligible losses, such states can be isolated through postselection of the 
SPDC output on the appropriate number of counts. In practice, the postselected events will 
also include other combinations of input photon numbers. Because typical detectors have 
limited or no photon number resolution, one needs to take into account also the deleterious 
contribution of higher total photon numbers. This leads us to the following definition of the 
parameter measuring the contamination of photon pairs with other terms that cannot be in 
general removed through postselection: 

.. = 1-,.^ (14) 

An analogous definition can be given for quadruples of photons, which ideally should be 
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FIG. 2: Contour plots of the contamination parameters for single photon pairs ei (upper plot) 
and double photon pairs £2 (lower plot) as a function of the overall efficiency r\ and the respective 
production rates pi^i and p2,2- The calculations have been carried out in the regime of single selected 
modes when M. = \. 



prepared in a state |2,2): 

£4 = 1- . (15) 

In Fig. [2] we depict contour plots of the contamination parameters £2 and £4 as a function of 
the production rates and the overall efficiency. It is clearly seen that the non-unit efficiency 
imposes severe bounds on the production rates that guarantee single or double photon pair 
events sufficiently free from spurious terms. The graphs also imply that strong pumping is 
not a sufficient condition to achieve high production rates, but it needs to be combined with 
a high efficiency of collecting and detecting photons. 
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FIG. 3: (Color online) Experimental setup. RegA, Regenerative amplifier; FL, IL, lenses; XSH, 
second harmonic crystal; DM, dichroic mirrors; BG, blue glass filter; ND, neutral density filter; 
HWP, half wave plate; X, downconversion crystal; RG, red glass filters; FC, fiber coupling stages; 
IF, interference filter; Dl, D2, avalanche photodiodes; CE, coincidence electronics. Thin solid lines 
indicate single mode fibers, dashed lines — multimode fibers. 



IV. EXPERIMENTAL RESULTS 

The experimental setup is depicted in the Fig. [3l The master laser (RegA 9000 from 
Coherent) produces a train of 165 fs FWHM long pulses at a 300 kHz repetition rate centered 
at the wavelength 774 nm, with 300 mW average power. The pulses are doubled in the second 
harmonic generator XSH based on a 1 mm thick beta-barium borate (BBO) crystal cut for 
a type-I process. Ultraviolet pulses produced this way have 1.3 nm bandwidth and 30 mW 
average power. They are filtered out of the fundamental using a pair of dichroic mirrors DM 
and a color glass filter BG (Schott BG39), and imaged using a 20 cm focal length lens IL 
on a downcoversion crystal X, where they form a spot measured to be 155 /im in diameter. 
The power of the ultraviolet pulses, prepared in the polarization perpendicular to the plane 
of the setup, is adjusted using neutral density filters ND and a motorized half waveplate. 
The type-I down-conversion process takes place in a 1 mm thick BBO crystal X cut at 29.7° 
to the optic axis, and oriented for the maximum source intensity. 

The down-converted light emerging at the angle of 2.5° to the pump beam is coupled 
into a pair of single mode fibers placed at the opposite ends of the down-conversion cone. 
The fibers and the coupling optics define the spatial modes in which the down-conversion 
is observed f3|. The coupled photons enter the loop detector {s] in which light from either 
arm can propagate towards one of the detectors through eight distinct paths. The minimal 
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delay difference between two patlis is 100 ns, more tfian twice tlie dead time of the detec- 
tors. Finally the photons exit the fiber circuit, go through interference filters IF and are 
coupled into multimode fibers which route them directly to single photon counting modules 
SPCM (PerkinElmer SPCM-AQR-14-FC) connected to fast coincidence counting electron- 
ics (custom-programmed Virtex4 protype board ML403 from Xilinx) detecting events in a 
proper temporal relation to the master laser pulses. The measurement series are carried out 
with pairs of interference filters of varying spectral widths. 

The measurement proceeds as follows: for each pair of interference filters the loop detector 
is calibrated using data collected at a very low pump light intensities, when the chance of 
more than one photon entering the fiber circuit is negligible compared to the rate single 
photons appear. This allows one to calculate the complete matrix of conditional probabilities 
Pk,n of observing k detector clicks with n initial photons j^, Q for one arm, and analogously 
P/^ for the second arm. The losses in the detectors and in the fiber circuit are assumed 
to contribute to the overall efficiencies r] and r]'. Thus Pk^n and P/„ describe lossless loop 
detectors for which Pi i = 1. After the calibration, the counts are collected for approximately 
10® master laser pulses for each chosen intensity of the pump, which yields the probabilities 
Pk,i of observing k clicks in one arm and I in the other one. These probabilities are related to 
the joint count probability /?„ ^ corrected for combinatorial inefficiencies of the loop detector 
through the formula: 



The probabilities pn,m can be retrieved from experimentally measured pk,i using the maxi- 
mum likelihood estimation technique jl^ . An exemplary joint photon number distribution 
reconstructed from the experimental data is shown in Fig. [H The results of the recon- 
struction are subsequently used to calculate the parameters of the source discussed in the 
preceding section. 

In Fig. [5] we depict the reconstructed equivalent number of modes A4 in single arm for 
different filtering and pump intensities. This quantity can be also understood as a measure 
of how many incoherent modes a photon from the source occupies. Naturally, drops 
with application of narrowband spectral filtering since it erases the information on the exact 
time the photon pair was born in the nonlinear crystal. It also is seen that Ai is practically 
independent on the pumping intensity, as predicted in Sec. Illli 




(16) 
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FIG. 4: An exemplary joint photon number distribution measured without spectral filters. The 
average photon numbers in the two arms are (n)=0.15 and (n')=0.18. 




FIG. 5: The equivalent number of modes M. in a single arm as a function of the average photon 
number (n), for measurements carried out without interference filters (circles), with 10 nm FWHM 
filters (squares) and 5 nm FWHM filters (crosses). 



The average overall efficiency of the squeezed state calculated as r/ = 1 — (5^) is shown in 
Fig. [6l It is seen that the efficiency decreases with an application of narrowband filtering, 
which again is easily understood. In addition, rj exhibits a very weak dependence on the 
pump intensity which again agrees with theoretical predictions for the regime when the 
average number of photons is much less than one. A relatively large difference in t] between 
10 nm and 5 nm interference filters can be explained by the fact that in the latter case the 
selected bandwidth becomes narrower than the characteristic scale of spectral correlations 
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FIG. 6: The equivalent detection efficiency r] as a function of the total average photon number 
(n), measured without interference filters (circles), with 10 nm FWHM filters (squares) and 5 nm 
FWHM filters (crosses). 

within a pair. Consequently, the filter in one arm selects only a fraction of photons conjugate 
to those that have passed through the filter placed in the second arm. This observation is 
consistent with the determination of the parameter Ai, which shows that for 5 nm filters 
effectively single spectral modes are selected. 

Finally, in Figs. [7] and [8] we plot the respective contamination parameters for single 62 
and double £4 photon pairs, that describes non-postselectable contributions of other photon 
terms to the output. As expected, the contamination grows with decreasing overall efficiency 
rj that corresponds to narrowing the spectral bandwidth, as well as with the increasing pair 
production rates. It is noteworthy that the measurement results agree well with fitted 
theoretical curves whose parameters match those that can be read out from Figs. [5] and [6l 
It is seen that in the regime the experiment was carried out the contamination is rather 
significant, and it would affect substantially effects such as two-photon interference. This 
limits in practice the energy of pump pulses — and consequently the production rates — that 
result in photon pairs sufficiently free from spurious terms. The effect is particularly dramatic 
in the case of £4, where genuine double pairs appear only in approximately half of all the 
events. 
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FIG. 7: (color online) The contamination parameter £2 for single photon pairs as a function of 
the pair production rate, measured without interference filters (circles), with lOnm FWHM filters 
(squares) and 5nm FWHM filters (crosses). The parameters of the fitted theoretical curves are: 
M = 16.6 and ry = 4.5% (red soHd line), M = 3.1 and ry = 4.2% (red dashed line) and M = 1.5 
and rj = 2.6% (red dash-dot line). 
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FIG. 8: (color online) The contamination coefficient £4 for double photon pairs as a function of 
the production rate, measured without interference filters (dots). The parameters of the fitted 
theoretical curve are: = 16 and 1] = 4.9% (red dashed line). Measurements with interference 
filters did not provide statistically significant data. 
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V. CONCLUSIONS 



Despite the simplicity of the basic concept, the application of SPDC sources in more 
complex quantum optics and quantum information processing experiments requires a care- 
ful choice of operating conditions. The actual output state is a result of a subtle interplay 
between the pump strength, spatial and spectral filtering of the output, and the losses expe- 
rienced by the signal. In this paper we concentrated on the features of the photon number 
distribution and developed an effective theoretical description. Starting from the fact that 
the SPDC output processed by arbitrary passive linear optics is given by a gaussian state, 
we introduced effective parameters characterizing the joint photon number distribution. The 
first parameter is the effective number of orthogonal modes that impinge on each of the de- 
tectors. This parameter carries information about the modal purity of the photons produced 
in each arm. The second parameter characterizes average losses experienced by the SPDC 
output, and can be calculated as a suitably normalized variance of the count difference. 
Finally, we also proposed and determined experimentally a measure of the contamination 
of photon pairs by other spurious terms that in general cannot be rejected by postselection 
and may contribute to an unwanted background. We showed that the presence of such a 
background puts stringent requirements on the detection efficiency if a bright source with 
high pair generation probability is desired. These constraints become very challenging when 
designing multiple-pair experiments. 
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